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Abstract. We report the field and temperature dependence of the low-temperature 
specific heat down to 400 mK and in magnetic fields up to 9 T of the electron-doped 
Ba(Feo.92Coo.o8)2As2 superconductor. Using the phonon specific heat obtained from 
pure BaFe2As2 we find the normal state Sommerfeld coefficient to be 18 mJ/mol K 2 and 
a condensation energy of 1.27 J/mol. The temperature dependence of the electronic 
specific heat clearly indicate the presence of the low-energy excitations in the system. 
The magnetic field variation of field-induced specific heat cannot be described by single 
clean s- or d-wave models. Rather, the data require an anisotropic gap scenario which 
may or may not have nodes. We discuss the implications of these results. 



PACS numbers: 71.20.Eh, 71.55.Ak, 72.15.Eb, 72.15.Jf, 75.50.Pp 



Submitted to: New J. Phys. 



Gap structure in the electron- doped Iron-Arsenide Superconductor 



2 



1. Introduction 

The recent discovery of superconductivity in Fe-based pnictides _RFeAsO[U [2] (i?-rare 
earth) has created a new era in superconductivity research and stimulated a great 
interest in these compounds [3j. Subsequently, other types of superconducting materials 
containing FeAs layers were discovered such as binary chalocgenides Fei +a; Se|H EJ, so 
called "111" compounds [61 [7] LiFeAs or NaFeAs and 122-systems AFe2As2 where A is an 
alkaline earth [HI [91 HO]. BaFe2As2, prototypical member of the latter family, crystalizes 
with the tetragonal ThCr 2 Si 2 -structure type and at ambient pressure exhibits structural 
and spin-density-wave (SDW) transitions at about 140 KpJ]. Suppression of the SDW 
state by either applied pressure [12] or chemical doping [8| [TO] results in superconductivity. 
Despite a large theoretical and experimental effort (see Ref.pSJ [H]) to understand the 
nature of the superconductivity in these materials there are still many open questions 
that have not been resolved such as the pairing mechanism and the symmetry of 
the order parameter. Moreover, the experimental results reported so far are often 
contradictory, ranging from nodal to fully gapped isotropic superconductivity. Even 
within the Co doped BaFe2As2 family the situation is unclear. While surface sensitive 
measurements such as ARPES[15] and STM[16] claim fairly isotropic gap values, 
penetration depth[T7], /xSR[18], NMRjJJ], thermal conductivity [20] EE] [22], specific 
heat|23j, and Raman scattering [21] argue that an anisotropic gap is necessary, although 
the details vary between these measurements as well. Often, the superconducting gap 
structure is discussed in terms of the s± model, with a sign reversal of the order 
parameter between different Fermi surface sheets [25 ] 126 ] [27] . However, the large sample, 
family, and doping dependence may favor scenarios where the low energy excitations, 
possibly nodal, depend strongly on the particular sample being studied and the probe 
used to investigate them (e.g. Ref . [28 | [29 | 130]). 

In this paper we present results of our detailed studies of the specific heat of the 
electron-doped Ba(Fe .92Co .o8)2As2 superconductor. By subtracting the lattice contri- 
bution (obtained from measurements on non-superconducting samples), we can extract 
the full electronic T-dependence. The temperature and magnetic field dependent data 
imply an anisotropic gap structure. 

2. Experimental details 

The large single crystals of Ba(Fe .92Coo.o8)2As 2 have been grown out of FeAs flux with 
the typical size of about 2x1.5x0.2 mm 3 . The samples crystalize as well-formed plates 
with the [001] direction perpendicular to the plane of the crystals. The doping level 
was determined by microprobe analysis. More details about the synthesis and char- 
acterization of the samples may be found in Ref. [9]. Based on electrical resistivity 
measurements, T c was established to be 20 K (zero resistance), in agreement with the 
heat capacity results presented here. The heat capacity was measured down to 400 mK 
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and in magnetic fields up to 9 T using a thermal relaxation method implemented in a 
Quantum Design PPMS-9 device. All data measured in field were field cooled. 

3. Results and discussion 

The temperature dependence of the specific heat of Ba(Fe .92Coo.o8)2As 2 is shown in 
Fig. 1. As can be seen from the figure, a pronounced anomaly of specific heat is observed 
at T c . A magnetic field of 9 T strongly suppresses the anomaly and moves it to lower 
temperatures. In general, the total specific heat of any system is the sum of several 
different excitations: 

C tot {T) = C el {T) + C ph (T) + C mag {T) + ... (1) 

where C e i(T), C p h(T) and C mag {T) describe electronic, lattice and magnetic 
contributions to the total specific heat, respectively. 




T[K] 

Figure 1. (Color online) The heat capacity of Ba(Feo.92Coo.os)2As2 measured in 
T (circles) and 9 T (squares). The solid line describes the normal state specific heat 
(see text). The inset shows the difference between measured C/T and the phonon 
contribution. The T c obtained by entropy balance construction is 20 K. 

In order to estimate the phonon contribution in our system, we assume that the 
phonon part of the specific heat is independent of doping. Thus, we determine the 
phonon specific heat from the parent compound BaFe2As2- BaFe2As2 exhibits a spin- 
density-wave transition at about 140 K so one may expect the presence of a magnetic 
contribution to the low-temperatures specific heat in the system. However, recent 
inelastic neutron scattering experiments show that, in the ordered state, spin-wave 
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excitations have a gap of about 10 meV (A « 116 K) [21] . Consequently, at temperatures 
below 40 K, C mag is negligible, and we separate the contributions to the specific heat of 
the parent compound as C = 7T + C p h- 

Thus, to describe the experimental data of Ba(Fe .9 2 Coo.o8)2As 2 above T c = 20 K, 
we adjust gamma to obtain the best agreement between the data and C = 7T + C p h- 
There are four points which give us confidence that our determination of 7 and C p h for 
the doped compound is reasonable, (i) the good agreement of C — 7T + C p h above 
20 K (see Fig.l). (ii) 7 which we obtain in this procedure provides accurate entropy 
balance for the electronic specific heat in the superconducting state below T c . (iii) 
the condensation energy which we obtain from the resulting analysis is quantitatively 
consistent with other measurements (see below), (iv) Calculations and inelastic x-ray 
scattering measurements indicate that the phonons below 10 meV are independent of 
dopping[3~2| 133] . 

The value of the normal state electronic specific heat may be compared with the 
density of states at the Fermi level calculated for pure BaFe2As2. Within a single band 
model approximation, 7 = 18 mJ/mol K 2 gives N(Ep) to be 7.64 eV _1 /f.u. Using LDA 
approximation together with GGA-PBE [31] or general potential in LAPW method [35J 
the calculated values of N(E F ), for the parent BaFe 2 As 2 , are 3.93 and 3.06 eV _1 /f.u., 
respectively. Thus the mass renormalization is roughly a factor of 2, consistent with 
mass renormalization determined by optics [36] and ARPES[37]. 

At 20 K the specific heat exhibits a jump AC/T C = 24 mJ/mol K 2 (see the inset 
in Fig.l) being consistent with previous reports (Ref. [9] [231 ESI EH])- It is also similar 
to AC/T C = 28 mJ/mol K 2 obtained for Ba . 6 Ko.4Fe 2 As 2 [23l HQ]. The size of the jump 
at T c depends on the details of the superconducting state. Taking 7 = 18 mJ/mol K 2 
the ratio AC /T c ^ = 1.33 is very close to, albeit smaller than, the weak-coupling BCS 
value 1.43. 

Using the normal state specific heat, we can extract the condensation energy 
and relate it to the thermodynamic critical field. This quantity can be obtained by 
integrating the entropy difference: 



where S n and S s denote entropy at the normal and superconducting state 
respectively. In the case of Ba(Fe .9 2 Co .o8)2As 2 this analysis give U = 1.27 J/mol and 
the thermodynamic critical field H c = 0.23 T. Using a penetration depth A = 325 nm 
obtained from MFMjJT] at a slightly different doping level and consistent with /xSR 
results [4"3] and coherence length £ = 27.6 obtained by STSfTB] gives a Ginzburg-Landau 
parameter K = j = 118. Using this value and our thermodynamic critical field H c we 
obtain H c2 = 38 T, from the expression H c2 = \/2KH c , in reasonable agreement with 
the published value of about 40 T|44j and the value of 39 T obtained from the slope of 
the upper critical field measured by specific heat and the expression if c2 =0.69^^T c [42]. 
This provides additional confidence in our phonon substraction 




(2) 
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Figure 2. (Color online) Temperature dependence of the electronic specific of 
Ba(Feo.92Coo.o8)2As2. The dashed and solid lines are theoretical curves based on BCS 
theory (Eq.3) with a one and two s-wave gaps, respectively (see text). 

We begin our investigation on the possible symmetry of the superconducting gap by 
examining the temperature dependence of the electronic specific heat. Fig. 2 displays the 
non-lattice part of the specific heat of Ba(Feo.92Coo.os)2As2 obtained by subtracting the 
phonon contribution together with a small Schottky contribution below 1 K (see below). 
At low temperatures, a sizeable residual specific heat coefficient 70 = 3.7 mJ/mol K 2 
is observed in this system. Similar behavior has also been reported in Bao.6Ko.4Fe2As2 
(70 = 7.7 mJ/mol K 2 )[l5] and Ba(Fei_ a Co x )2As2 (70 ~ 3 mJ/mol K 2 for optimal doped 
samples) [23] . Interestingly, the sizeable value of the residual specific heat coefficient 
also has been observed in superconducting cuprates[l6l 07] . For Ba(Feo.92Co .o8)2As 2 , 
70 = 3.7 mJ/mol K 2 amounts to 20 % of 7. We rule out that 70 originates from 
a non-superconducting portion of the sample based on the fact that x-ray analysis 
limits impurity phases to less than 5%, and we observe full diamagnetic shielding from 
magnetization measurements. Alternative explanations for the origin of the residual 70 
include pair breaking effects of an unconventional superconductor [56] . crystallographic 
defects or spin glass behavior. 

To fit the C e i/T data in figure 2 we use the BCS expression for specific heat: 

Cbcs = tj t J™ dy (-^) [flnf + (1 - /) In (1 - /)] (3) 
where t = f is the Fermi function / = — ^ — , E = y/e 2 + A 2 and V = \ 

e k B T +1 

(see Ref.[l8]). 
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Figure 3. (Color online) The heat capacity of Ba(Feo.92Coo.os)2As2 plotted as CJT 
vs T 2 with magnetic field applied along the c-axis. 

The results of our analysis using a single s-wave gap|49j (dashed blue line) and 
two separate s-wave gaps (solid red line) are shown in Fig. 2. In an s-wave model a 
residual linear term must be an extrinsic contribution to the specific heat, and hence 
we subtracted 70 from 7 for the purposes of these fits. The normal state Sommerfeld 
coefficient 7 = (18 - 3.7) mJ/mol K 2 and T c = 20 K were held fixed during the fits. 
The gap value obtained from the single gap fit is 3 meV and may be compared with 
A = 6 meV derived for hole-doped Bao.6Ko.4Fe2As2 with T c « 37 Kjl5]. Taking 
A = 3 meV and T c = 20 K gives A /(/csT c ) = 1.74, close to the weak coupling 
value of 1.76. However, as can be seen from Fig. 2 the single gap fit does not describe 
the data sufficiently and clearly indicates the presence of low energy excitations in 
the system below 8 K. A much better description is obtained by fitting the data to 
C = (I-A)Cbc's(Ai) + ACbcs(A 2 ) which gives the solid red line and the parameters 
Ai = 1.65 meV, A 2 = 3.75 meV and A = 0.62. While this fit provides a reasonable 
description of the data, we emphasize that there are multiple anistropic gap descriptions 
that could provide a similarly good fit. From this data alone we cannot determine 
whether or not nodes exist. 

The low-temperature specific heat of Ba(Feo.92Coo.os)2As2 measured in several 
magnetic fields is presented in Fig.3. Below 1 K an upturn in C/T is observed. Such 
behavior has been already observed in Ba(Feo.92Coo.os)2As2 samples grown by In flux as 
well as in Ba .6K 4 Fe 2 As 2 and FeSe single crystals[50ll4"5"|l51]. The origin of the anomaly 
is most probably related to the presence of a small amount of magnetic impurities. With 
increasing magnetic field the upturn shifts to higher temperatures and transforms into 
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Figure 4. (Color online) Field-induced change in low temperature specific heat of 
Ba(Feo.92Coo.o8)2As 2 , obtained at K by extrapolating the experimental data of Fig. 3 
to zero temperatures (see text). The green dotted line and red dashed line represent 
field dependencies expected for s-wave and rf-wave descriptions, respectively. The blue 
solid line is a theoretical curve for anisotropic s-wave superconductors (see text). 

a maximum, which gradually broadens and diminishes in magnitude with further rising 
field. Such behavior indeed resembles that expected for a degenerate ground state split 
due to the Zeeman effect in internal and external magnetic fields. Thus, to avoid the 
effects of this magnetic contribution, we rely on a linear extrapolation of the data from 
above 1.5 K. The red solid lines in Fig. 3 display the linear tendency of low-temperature 
specific heat displayed as C/T(T 2 ). We have used it to determine the electronic specific 
at K by extrapolating the low-temperature data to zero temperature. 

The so obtained A7 = ^- c ^~ c ^ derived at K, as described above, as well as at 
several other temperatures is presented in Fig. 4. Similar data were obtained in Ref. [23J. 
The data have been presented in the form A7/7 vs H/H c2 with 7 = 18 mJ/mol K 2 and 
H C 2 = 39 T. In fully gapped superconductors the localized quasiparticle states in vortex 
cores result in A7 proportional to H/H C 2 since the number of vortices is proportional 
to the magnetic field. As can be seen from the figure, the specific heat rises much faster 
with field than expected for a simple s-wave gap. This could indicate an anisotropic gap 
or point to a field dependent coherence length [52| [53] . The latter scenario, however, is 
unlikely to be a sole explanation as the temperature dependence indicates low-energy 
excitations in the system (see Fig. 2). Thus, we further explore the expectation of an 
anisotropic gap for an explanation of the A'y(H) dependence. 
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It was shown theoretically [55J by using microscopic quasiclassical theory that 
an anisotropic gap structure can display a significant field dependence in A7(H) . 
Taking into account this model the experimental data may be well described at low 
magnetic field (see the blue solid line in Fig. 4) using the gap anisotropy ratio a = 0.5 
(A m j n / A max = 0.5), which is in resonable agreement with the two gap fit of the 
temperature dependence in Fig. 2. 




Figure 5. (Color online) Normalized field-induced change in low temperature specific 
heat of Ba(Feo.92Coo.o8)2As2, obtained at K, versus magnetic field. The red solid line 

is a fit of Eq.4 to the experimental data. Inset: the same data plotted versus a (^jt~^ 
with a = 0.5 and H c2 = 40 T. 

An extreme anisotropic gap is that of a nodal superconductor as occurs in the 
cuprates. For clean d-wave superconductors, it has been shown by Volovik[54"] that 
the quasiparticle excitation spectrum is shifted by the superfluid velocity, resulting in 
A7 oc ^/H/H c2 . The quasiparticles that contribute to the density of states are coming 
from regions far from the vortex cores and close to the nodes. This situation is presented 
in Fig. 4 by the red dashed line. As can be seen from the figure A7 is not increasing so 
strong as expected for clean d-wave superconductors. 

Finally, we consider an additional alternative. The finite density of states 70 
observed in Ba(Feo.92Coo.o8)2As 2 may also be related to impurity scattering effect 
as expected for a dirty d-wave superconductor. It has been shown by Kiibert and 
Hirschfeld|56] that in the dirty d-wave limit, A7 behaves like HlogH at the lowest fields 
{H c i < H <C H C 2). In this approach the field dependence of field-induced specific 
heat may be expressed by [57]: 

A 7 A (H\ 1 _JB\ 
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where A = 0.322 {Af) 2 , B = ^^f, a w 0.5 and Ao and T denote a maximum 
of the gap and impurity scattering rate, respectively. The solid line in Fig. 5 is a fit 
of Eq.4 to the experimental data, resulting in the parameter ^ = 78. During this 
analysis the value of H C 2 has been fixed to 40 T. The inset in Fig.5 shows the data for 
Ba(Fe .92Coo.o8)2As 2 together with a curve calculated for a clean d-wave superconduc- 
tor. As may be seen the field dependence of r y(H)/ r y n as well as the magnitude of the 
residual specific heat may be consistent with the dirty ci-wave scenario. The point of 
these two fits (Fig.4 and 5) is to demonstrate that we can not distinguish between these 
two different interpretations. The quality of agreement between the two gap analysis in 
Fig.4 and the dirty c?-wave analysis of Fig.5 is comparable. 

From the above analysis, we can determine that a single isotropic s-wave gap is in- 
capable of describing the specific heat data. The form of the anisotropic gap, however, 
cannot be determined by our data alone. We see that two extreme cases of multi- 
band s-wave, and a dirty d-wave scenario are each in reasonable agreement with the 
data. Results from NMR[T9], thermal conductivity PS ED E2], ARPESpS], penetration 
depth [17], /^SR[18] and Raman [24] also conclude that an anisotropic gap is necessary 
to describe their data on Co-doped BaFe2As2, although the extent to which varies from 
gapless to mild multiband behavior. A nodal gap imposed by symmetry, as in the case 
of the cuprates, is ruled out by the vanishingly small residual linear term of the thermal 
conductivity. Hence, the dirty rf-wave analysis applied above should not be directly 
applicable. However, accidental nodes as anticipated in some spin-fluctuation models 
of the pnictides (see Ref.[28| |29| 130]). cannot be ruled out by the thermal conductivity 
results [58] and are also consistent with our specific heat data. The accidental node sce- 
nario has the favorable aspect that the nodes could be lifted by disorder and/or doping 
which would help reconcile some of the seemingly contradictory results [59]. Further 
measurements as a function of doping and disorder are necessary to help elucidate the 
gap structure, not to mention the variations between different families and dopant atoms. 

4. Summary 

In summary, we have used low-temperature specific heat and its magnetic field re- 
sponse to explore details of the symmetry of the superconducting gap in electron-doped 
Ba(Fe .92Co .o8)2As2 superconductor with T c = 20 K. Using the phonon part of the 
specific heat of pure BaFe2As2, we determine the normal state Sommerfeld coefficient 
in Ba(Feo.9 2 Coo.o8)2As2 to be 18 mJ/mol K 2 . The temperature variation of the elec- 
tronic specific heat below T c may be well described by the presence of two supercon- 
ducting gaps, pointing to complex gap structure in the system. The field-induced low- 
temperature specific heat can not be explained by simple clean s- or c?-wave descriptions. 
Its behavior also indicates a strongly anisotropic gapped superconductor. 
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Note: During completion of this manuscript we became aware of ref. [SD] which 
used a similar procedure to determine the phonon contribution to the specific heat of a 
similar crystal. The resulting temperature dependence was analyzed within a two band 
model with results in good agreement with ours. 
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